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CN ■ ABSTRACT 
>• ■ 

The physical processes involved in diffusion of Galactic cosmic rays in the interstellar medium are 
addressed. We study the possibility that the nonlinear MHD cascade sets the power-law spectrum of 
\ turbulence which scatters charged energetic particles. We find that the dissipation of waves due to 

■ the resonant interaction with cosmic ray particles may terminate the Kraichnan-type cascade below 
^—1 ' wavelengths 10^'^ cm. The effect of this wave dissipation has been incorporated in the GALPROP 

, numerical propagation code in order to asses the impact on measurable astrophysical data. The 

■ energy-dependence of the cosmic-ray diffusion coefficient found in the resulting self-consistent model 
' may explain the peaks in the secondary to primary nuclei ratios observed at about 1 GeV/nuclcon. 

Subject headings: diffusion — MHD — elementary particles — turbulence — waves — cosmic rays 

o ■ 

H ■ 

+-? ■ 1. INTRODUCTION 

C/2 . 

■ The Galactic cosmic rays - the gas of relativistic charged particles with high energy density - cannot always be 
^ '■ treated as test particles moving in given magnetic fields. In particular, the stochastic acceleration of cosmic rays by 

• I— I MHD waves is accompanied by the damping of the waves, since the wave energy is dissipa ted. The rate of wave 
' damping on cosmic rays through the cyclotron resonance interaction was first estimated bv iTidnianI (jlQBfif l. If we 
?H \ exclude cold HI regions (where the waves are damped by collisions of ions with neutral atoms) and also perhaps 
. 5t I regions of the interstellar gas with very high temperature T ~ 10^ K and weak magnetic field (where Landau damping 
on thermal particles is high), then one finds that this mechanism of dissipation could dominate over other known 
mechanisms in the interva l of wavelength s 10^^ cm to 10^^ cm. The non-resonant interaction of diffusing cosmic rays 
with magnetosonic waves lIPtuskinlUQSl is important only at large wavelengths and is not important for the present 
investigation. We also do not consider the large body of instabilities in cosmic rays that may arise because of the 
non-equilibrium distribution of charged energetic particles caused by possible strong anisotropy or large gradients and 
which may amplify waves in the background plasma - see e.g. iBer ezinskii ct al. (1990) and Diehl ct al. (2001) for 
a review of such processes. Cyclotron wave damping on cosmic rays changes the wave spectrum in the interstellar 
medium, which in turn affects the particle transport since the cosmic ray diffusion coefficient is determined by the 
level of turbulence which scatters charged particles. Thus in principle the study of cosmic ray diffusion requires a 
self-consistent approach. We shall see below that the effect of cosmic rays on interstellar turbulence should be taken 
into account at energies below a few GeV/nucleon, and this may result in a considerable increase in the particle 
diffusion coefhcient. This picture will be also verified by its consistency with observations of interstellar turbulence. 
Consistency of this picture with observations of interstellar turbulence is also discussed 

It is remarkable that the interpretation of cosmic ray data on secondary nuclei may require this effect. The secondary 
nuclei are produced in cosmic rays in the course of diffusion and nuclear interactions of primary nuclei with interstellar 
gas. The ^H, '^He, Li, Be, B and many other isotopes and elements are almost pure secondaries. Their abundance 
in cosmic ray sources is negligible and they result from the fragmentation of heavier nuclei. Cosmic ray antiprotons 
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and the major fraction of positrons are among the secondary species as well. The ratio of secondary to primary 
nuclei such as B /C has a peak at about 1 GeV/n and decreases both with increasing and decreasing energy, e.g. 
lEiigel mann "eTall lfT990|) . The high energy behavior is naturally explained by particle scattering on interstellar MHD 
turbulence with a power-law spectrum in wavelength, but the required sharp increase of the diffusion coefficient at 
small energies has usually been considered as improbable since it needs a drastic "p hysically unjustified" bending down 
of the wave spectrum at wave numbers k > 1/3 x 10"'^^ cm~^, see e.g. I.Tones et al.. ()2001.) . .Moskalenko et al. (2002). 
Alternative explanations of the peaks in secondary to primary ratios without invoking peculiarities in the cosmic ray 
diffusion coefficient have been suggested. The two most popula r of these involve diffusive-convective particle transport 
in the hypothetical Galact ic wind ( Jones Il979t IPtu skin et cd.|[r997) and the stochastic reacceleration of cosmic rays 
by interstellar turbulence l) Simon et allilQSfit ISeo fc Ptusk iTri994'l. These processes are incl uded in the numerical 
computations of cosmic raj^ransport in the Galaxy in the frameworks of the GALPROP code ijStrong fc M oskalenkol 
119981 ^oskalenko et al. 2002) . The situation with the interpretation of the energy dependence of secondaries in cosmic 
rays is still uncertain. The peaks in the secondary to primary nuclei ratios calculated in diffusion-convection models 
seem to be too wide. Models with reacceleration re produce the shape of the peaks but the absolute flux of antiprotons 
turns out to be too low compared to observations iMos kalenko et al.ll20 0^. This is why the effect considered in the 
present work is so important. We investigate a self-consistent model of cosmic-ray diffusion in interstellar turbulence 
where the wave damping on energetic particles is taken into account, and implement this effect in the GALPROP code. 
To make the first calculations reasonably tractable we keep only the most essential features of cosmic-ray diffusion in 
random magnetic fields and use the simplest description of the nonlinear wave cascade in the interstellar turbulence. 

2. EQUATIONS FOR COSMIC RAYS 

The steady state transport equation that describes d iffusion and convective transport of cosmic ray protons and 
nuclei in the interstellar medium is of the form (see e.g.. lBerezinskii et alJll990l for discussion): 



d 

Vi:'V* + V(u*) - — 
op 
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Here 5'(p, r) is the particle distribution function in momentum p normalized on total cosmic ray number density as 
A'cr = / dp^, D{p,r) is the spatial diffusion coefhcient, K{p,r) is the diffusion coefficient in momentum, u(r) is the 
velocity of large-scale motions of the interstellar medium (e.g. the velocity of a Galactic wind), pioss — dp/dt < is 
the momentum loss rate for the energetic charged particles moving through the interstellar medium, T(r) is the time 
scale for nuclear fragmentation. If needed, the supplementary term which describes radioactive decay can be added to 
eq. The spatial boundary condition for is ^' |s= 0, corresponding to the free exit of cosmic rays from the Galaxy 
to intergalactic space where their density is negligible. The region of cosmic ray diffusion is a cylinder of radius 30 kpc 
and total width 2H (H w 4 kpc). The source term q{p,v) includes both the direct production of primary energetic 
particles accelerated from the thermal background in Galactic sources (e.g. supernova remnants) and the contribution 
to the nuclei considered via the processes of nuclear fragmentation and radioactive decay of heavier nuclei. The typical 
value of the diffusion coefficient found from the fit to cosmic ray and radioastronomical data is D ~ 3 x 10^* cm^/s at 
energy '^l GeV/n, giving a diffusion mean free path I = ZD /v ~ 1 pc (w ~ c is the particle velocity). 

On the "microscopic level" the spatial and momentum diffusion of cosmic rays results from the particle scattering 
on random MHD waves and discontinuities. In the linear approximation, the Alfven (with the dispersion relation 
Ci;(k) = the fast magnetosonic (w(k) = and the slow magnetosonic (w(k) ~ fc||V^) waves can propagate in 

a low f3 plasma, /3 = (14/Ki)^ < 1, where 14 — B / ^Airp is the Alfven velocity and Vg = \J P^j p is the collisionless sound 
velocity (for MHD conditions it is a factor of a square root of the gas adiabatic index larger), where B, Pg and p are the 
magnetic field strength, the pressure, and the mass density of the interstellar gas respectively. In addition to waves, 
the static entropy variations can exist in the interstellar medium. In a collisionless plasma, the slow magnetosonic 
waves are not heavily damped only if the plasma is non-isothermal and the electron temperature considerably exceeds 
the ion temperature. The Alfven velocity is 1.8x 10^-B^g/v^ "^i^/Si where n is the number density of hydrogen atoms. 
The value of the sound velocity is 7.7 x Wf'^/Tl cm/s in neutral gas with temperature T = 10^24 K; is larger in 
the ionized gas because of the free electron contribution to gas pressure Pg. The approximation of low /3 plasma is 
valid in the dominant part of the interstellar medium at P w 5 pQ since typically n = 0.002 cm^'^ and T — 10^ K in 
hot HII regions, n = 0.2 cm""^ and T = 8 x 10'^ K in warm intercloud gas, n = 30 cm~^ and T = 100 K in clouds of 
atomic hydrogen, n = 200 cm~'^ and T = 10 K in molecular clouds and thus /3 « 0.1 — 0.3 everywhere. The effective 
"collision integral" for energetic charged particles moving in small a mplitude random field s SB < S! P can be taken 
from the standard quasi- linear theory of plasma turbulence, see e.g. iKennel fc EngelmannI l|196fif l. The "coUisional 
integral" averaged over fast gyro-rotation of particles about the magnetic field B contains the spectral densities of the 
effective frequencies of particle collisions with plasmons in the form 



^^"(k,s,p) = -Tyn^ [(1 - '5«o)5(fc|| - s/(rg |M|)M^(k)(j2+i + Jl_^) 



+ _ fciv;/(z; ImD) Af|f (k)fcf fcl^]^ (2) 

fell > (as presented bv lPtuskinlll989(l . Here the waves are assumed symmetric about the magnetic field direction and 
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to have zero average helicity. The index a characterizes the type of waves including the direction of their propagation, 
/i is the cosine of pitch angle, the energy densities of random magnetic fields for perpendicular and parallel to the 
average field B components are Af"(k) and Mj"(k) respectively, Jm = Jm{k±rg^Jl — /i^) is the Bessel function. It 
is explicitly taken into account that the scattering at s = occurs only on the fast magnetosonic waves propagating 
almost perpendicular to the magnetic field. The Larmor radius is — pc/{ZeB) — 3.3 x 10^^i?Gv/-S^G cm, where 
the particle magnetic rigidity R = pc/Ze is measured in GV, and the average magnetic field is measured in ^G. 
The wave-particle interaction is of resonant character so that an energetic particle is predominantly scattered by the 
irregularities of magnetic field SB that have the projection of wave vector on the magnetic field direction equal to 
fc|| = ±s/ (rgii). The integers s = 0, 1,2... correspond to the cyclotron resonances of different orders. The efficiency 
of particle scattering depends on the polarization of the waves and on their distribution in k-space. The first-order 
resonance s = 1 is the most important for the isotropic and also for the one-dimensional distribution of random 
MHD waves along the average magnetic field, k || B. In some cases - for calculation of scattering at small ^ and for 
calculation of perpendicular diffusion - the broadening of resonances and magnetic mirroring effects should be taken 
into account. 



The evolution of the particle distribution function on time-scales At ^ 



and distances Az vv can be 



described in the diffusion approximation with the following expressions for the spatial diffusion coefficient along the 
magnetic field and the diffusion coefficient in momentum Dpp: 



(3) 



Dpp{p) 
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^ (4) 
see lBerezinskii etall ()199C1|1 where the analogous equations were derived for one-dimensional turbulence with k± = 0. 
Here Vf"(k) = a;"(k)/fc||, the summations contain integrals over k-space, and the terms which contain f^"(k, s — 0,p) 

should be corrected in eqs. lO, Q compared to eq. multiplied by (1 — ^^)^ in (E^^?) multiplied by (1 — /^^) 

in other terms. One can check that the interaction of energetic particles with slow magnetosonic waves is relatively 
weak (as /V^ <^ 1) and can be ignored. 

Locally, the cosmic ray diffusion is anisotropic and occurs along the local magnetic field because the particles are 
strongly magnetized, rg <^ I. The isotropization is accounted for by the presence of strong large-scale (~100 pc) 
fluctuations of the Galactic magnetic field. The problem is not trivial even in the case of relatively weak large scale 
random fields, since the field is almost static and the strictl y one-dimensional diffusion a long the magnetic fie l d line s 
does not lead to non-zero diffusion perpendicular to B, see iChuvilgin fc PtuskinI l)1993D . iGiacolone fc Jokipiil l)1999D . 
ICassc ct al. ( 2001). 

The eqs. Q and Q are too cumbersome for our p r esent appl ication. Based on eg. 1311 and with reference to the 
detailed treatment of cosmic ray diffusion bv lToptvginl ||198!tD and lBerezinskii et alJ l[T99flil . we use below the following 
simplified equation for the diffusion coefficient: 



D = vrgSV [127rfcrcsVK(fcrcs)] 



(5) 



1/rg is the resonant wave number, and W{k) is the spectral energy density of waves normalized as 

res <C B. 



where fc^os 

/ dkW{k) = SB'^/Att. The random field at the resonance scale is assumed to be weak, SB, 
The cosmic ray diffusion coefficient is equal to 



D 



B' 



3(1 - a)kl: 



SBl 



(6) 



for particles with rg < k^^ under the assumption of a power law spectrum of turbulence W{k) oc l/fc^~°, k > k^. We 
introduced here the principle wave number of the turbulence k^ and the amplitude of random field SB^ at this scale. 
The diffusion has scaling D cx v{j)/Z)°-. 
The diffusion in momentum is described by the following equation which is a simple approximation of eq. 

K^p'V^liQD). (7) 

Eqs. (0-0 give estimates of particle diffusion in position and momentum needed in eq. 10. They reflect the most 
essential features of cosmic-ray transport: the frequency of particle scattering on random magnetic field is determined 
by the energy density of this field at the resonance wave number /cres! the acceleration is produced by the waves 
moving with typical velocity T4 and is stochastic in nature. Eq. 10 implies equal intensities of waves moving along the 
magnetic field in opposite direction s, the imbalanced part of the total wave energy density should not be taken into 
account when calculating K{p), see Berezinskii et al.l ()1990() for details. Eqs. ©-0 are also valid for small-amplitude 
nonlinear waves including weak shocks. The equations can be used for the isotropic distribution of Alfven and fast 
magnetosonic waves and they give correct order of magnitude estimates for the wave distribution concentrated around 
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the direction of average magnetic field. It should be pointed out however that the isotropization of the diffusion tensor 
does not occur in the case of a pure parallel propagation of waves (k || B). Another special case is 2D turbulence 
with perpendicular propagation of waves (k _L B). In this case, the scattering occurs only for magnetosonic waves 
through the resonance s = which leads to a very large diffusion coefhcient, about a factor (u/T4)^ larger than given 
by eqs. (0), ©• 

Eq. ^ shows that the level of interstellar turbulence that is needed to account for the diffusion of GeV cosmic rays 
is very small: SB^^^/B^ = (1 — a)rg/l ^ 10^^ at k~J. ~ lO^'^ cm (if a ^ 0.5). An extension to smaller wave numbers 
gives (5B2(> k)/B'^ - lO-^{10^'^k)^-'' where k is in cm"!. 

3. EQUATIONS FOR INTERSTELLAR TURBULENCE 

The description of MHD turbulence is a complicated and not completely solved problem ev en in the case of small- 
amplitude r andom fields. Comprehens i ve re views of MHD turbulence have been given bv iVernial l)2004f) . and by 
iZhou et ai] l|2004H . lEhnergreen fc Scald 1)20041) . and lScalo fc ElmergreenI l)2004|) with apphcation to interstellar turbu- 
lence. 

The classic problem is the determination of the wave spectrum in the presence of sources at small wave numbers 
fe ~ /cl and the strong dissipation at much larger wave numbers (in some cases the cascade is inverse). Note that 
the spectrum of interstellar MHD turbulence determines the transp ort coefficients eqs. 0. According to the 
Kolmogorov-Obukhov hypothesis l)Kolmogorovlll94'll lObukhovlll941|) . the resulting spectrum at intermediate k, i. e. 
in the inertial range, is characterized by a constant energy flux to higher wave numbers. The hypothesis was origi- 
nally suggested for the description of developed hydrodynamic turbulence in incompressible fluids. The Kolmogorov 
spectrum is of the form W (k) cc The spectrum of weak acoustic turbu lence W(k) oc fc^^/^ was found by 

iZakharov fc SagdeevI ()197nD . This result was criticized by Kadom tsev fc Petvias hvili (1973). They argued that the 
developed shocks produce an additional dissipation of acoustic wave energy. The shock-dominated turbulence (also 
called the Burgers turbulence) is characterized by the spectrum k~^ irrespective of the nature of dissipation at the 
shock front. 

The presence of magnetic field in MHD turbulence complicates the issue because the turbulence becomes anisotropic 
and new types of waves arise in a magnetized medium. Generally, all gradients are larger perpen dicular to the field and 
all perturbations are elongated along the magnetic field direction even with isotropic excitation. Ilroshnikovl l) 196.11) and 
iKraichnanI l)1965(l gave the first phenomenological theory of MHD turbulence and obtained the spectrum W{k) cc k~^/^. 
They assumed that small-scale fluctuations are isotropic, which is contradictory because wave interactions break the 
isotropy in the presence of an external magnetic field. 

Direct observations of MHD turbulence in the solar wind plasma where (5^1 have shown the existence of a 
Kolmogorov-type turb ulence spectrum that contains waves moving both along the magnetic field and in almost per- 
pendicular directions l)Saur Bieberlll999D. Such a spectrum w as obtained over several decades of wave numbers 
in solar wind radio propagation studies iWoo fc A rmstrong 1979'). Numerical simulations of incompressible {fi ^ 1) 
MHD turbulence favored the Kolmogorov spectrum (Vcrma ct al, 1996) . 

Over the past decade, there has been a renewed interest in u nderstanding of magnetohydrodynamic turbulence as 
it apphes to interstellar magnetic field and density flu ctuations llGoldreic h fc Sridhar"199f; . 1997 : lNg fc Bhattacharied 



^97; Galt ier et al.l2000l see also earlier work by Shcba hn et alJl983j) . iGoldrcich fc Sridhar ( 199E ) exploited anisotropy 
in MHD turbulence and obtained Kolmogorov-like spectrum for the energy density of Alfven waves. The "elementary 
interactions" between Alfven waves satisfies the three-wave resonance conditions. However there is no exact relation 
between wave number and frequency in this case of strong turbulence. The main part of the energy density in this 
turbulence is concentrated perpendicular to the local magnetic field wave vectors fc_L « fc, while the parallel wave 

1/2 

numbers are small: fc|| ~ \kW{k)/ (i?Q/47r)] k±. The cascade is anisotropic with energy confined within the cone 

fc|| a k^/^^ . Numerical simulations have confirmed this concept (|Cho fc Vishniac|l2000|) . 

Although the formalism has been developed for incompressible MHD turbulence . iLithwick fc G oldreich' ("2001') argued 
that the compressibility does not essentially alter the results on the Alfven wave spectrum. The distribution of slow 
magnetosonic waves passively follows that of Alfven waves. The fast magnetosonic waves have an independent nonlinear 
cascade which is isotropic and has a Kraichnan-type spectrum W{k) oc fc""^/^. These conclusions were supported by 
numerical simulations by Cho fc Lazarian (2002). 

The description of weak MHD turbulence in low /3 plasma outlined above is probably not complete and needs further 
analysis before it is accepted as a standard model of interstellar turbulence. First, there is still the discrepancy between 
theoretical results of different authors. Thus co nsidering the scattering of Alfven waves and fast magnetosonic waves 
on slow magnetosonic waves, iKuz netsovl ()200lD found the Kraichnan-type spectra for all these types of waves with 
their preferentially parallel propagation, which disagrees with Lithwick fc Goldreich (2001J. Second, consideration of 
processes in turbulent coUis ionless plasmas at the kinetic level involves additional no nlinear processes of induc ed wave 
scattering on thermal ions (jLivshits fc Tsvtovich 1970') that may change the spectra (.Chas hei fc Shi shovl 19851) . Third, 
the real turbulence can be strongly intermittent, imbalanced, etc. ('e.g. ILithwick fc G^ldreid]ll2003ir which mav also 
affect the interstellar MHD spectrum. 

Information on the extended interstellar turbulence spectrum has been obtained from radio scintillation and re- 
fraction observations (sensitive to fluctuations of thermal electron density), measurements of the differential Faraday 
rotation angles from distant sources (mainly produced by fluctuations in the interstellar magnetic field), and the ob- 
servations of random motions in the interstellar gas. These data are consistent with the assumption that a single 
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close-to-Kolmogorov spectrum extends from scales 10^ to 10^° cm (|Lee &: Jokipiil[T976(l . see lArmstrong et al.l ()1981l 
11995ft and references therein. 

In the absence of an easily manageable and commonly accepted exact equation for the energy density W{k), we 
employ below the simplest steady state phenomenological equation that represents the concept of a wave cascade in 
the inertial range of wave numbers: 

d_ / kW{k) \ 

dk \ r„i ) 

ljTij|ll988t INorman fc Ferrarall 19961 this approach goes back to IChandrasekhailll948l and IHeisenberelll948|) . Here the 
approximation of a characteristic time T^i is used to account for the non-linear wave interactions which provide the 
transfer of energy in k-space**. 



-0 (8) 



Formally, the Kolmogorov-type spectrum W{k) cc k ^/■^ follows from eq. |(SJ) if Tni = TX = CAky^kW{k)/ (47r/ci) , 
where the constant Ca ~ 0.3 as can be estimated from the simulations of Verma et al (1996). The Kraichnan spectrum 
W{k) cc is obtained from eq. ^ if Tni = Tm = [Cuk'^Wik)/ (pVa)] "\ where Cm ^ 1. 

The Kolmogorov-type nonlinear rate 1 /T\ is relatively high since it is proportional to the amplitude of weak ran- 
dom field whereas the Kraichnan rate 1/7m i s proportional to the amplitude squared. Our preliminary estimates 
(jMoskalenko et al]l2003at iPtuskin et aUlWoSal) showed that the MHD cascade with the Kolmogorov rate is not sig- 
nificantly affected by damping on cosmic rays even if the waves propaga te along the magnetic field which makes the 
wave-particle interactions the most efficient. In addition, if the concept of lGoldrei ch fc Sridhar (1995') works for inter- 
stellar turbulence, then the Kolmogorov rate refers to the Alfven waves which are distributed almost perpendicular to 
an external magnetic field and thus do not produce any significant scattering of cosmic rays, see Section 2. So, we do 
not consider the damping of cascades with the Kolmogorov rate in what follows. 

For a cascade wit h the Kraichnan rate 1/Tm which describes an isotropic or close to parallel distribution of MHD 
waves (according to lGoldreich fc Sridhail ll995. it refers only to fast magnetosonic waves but not to Alfven waves), the 
equation for wave energy density can be written as follows: 



o, > ,-k^W'{k)\=-2Tik)Wik) + SSik-ki^), (9) 
dk / 

k > k]^. Here F = Fcr + Fth is the wave attenuation rate on cosmic-ray particles (Fcr) and on thermal particles 
(Fth), S characterizes the source strength. The main sources of interstellar turbulence are supernova explosions, winds 
of massive stars, superbubbles, and differential rotation of the Galactic disk. These sources produce strong random 
magnetic fields and probably initiate nonlinear wave cascades at the scale ^ ~ 100 pc. 

Below we ignore the dissipation on thermal particles (see discussion below) and set Fth = to study purely the effect 
of damping on cosmic rays. We actually deal with a small part of the possible global wave spectrum at wave numbers 
k > 10"^'' cm~^ which resonantly scatter cosmic rays with energies less than about 100 GeV/n. The existence of a 
single interstellar MHD spectrum over 12 orders of magnitude is an unsolved problem in itself and is beyond the scope 
of the present work. The quantities kj^ and W{ki^) are used here o nly for purposes of nor malization. 

The equation for wave amplitude attenuation on cosmic rays is ijBerezinskii et a'l]ll990|) : 

where pres (k) = ZeB/ck. 
The solution of eqs. lO, (|10|l allows us to find the wave spectrum: 



W{k) 



-3/2 



p™(fci) 



(11) 



k > /cl- Here the spectral wave density at the principal scale is determined by the source strength: W{ki^) = 

\/ pVs.S/Cm k^,^^"^ ■ The second term in square brackets is increasing with k. The wave damping on cosmic rays 
decreases the wave energy density and can even terminate the cascade if the expression in square brackets reduces to 
zero at some fc = fc*. The wave energy density should be set to zero at fc > fc* in this case. 
Now with the use of eqs. lO, (|ll|l one can determine the cosmic ray diffusion coefficient, which is 



D{P) = . 1/2 ■ (12) 



2CMB^r„v p^r.^.h'-^l'^ 



'WpcC^r^'Wik^) 



^ See also ICranmer van Ballegooiiei] 1200.?!^ and references therein for possible modifications of eq. j^J when a diffusive flux of wave 
energy in fc-space is assumed and when the equation is written for a power spectrum VF(k) not averaged over direction (here rS^kWi}^) = 
dkWik)). 
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p < PL, PL = Pros(fcL)- Here Do{p) = vrg^^ B'^ /[12Trk'[^^W{ki,)] cc is the diffusion coefficient calculated for 

the Kraichnan-type spectrum without considering wave damping. The second term in the denominator of eq. H12|l 
describes the modification of the diffusion coefficient due to wave damping. The effect becomes stronger to smaller p. 
The diffusion coefficient should be formally set to infinity at p < if the square bracket in eq. goes to zero at 
some p = p^,. The constant g characterizes the strength of the effect for a given cosmic ray spectrum ^'(p). 

As the most abundant species, the cosmic ray protons mainly determine the wave dissipation. Thus, with good 
precision only the proton component with Z = 1 needs to be taken into account to calculate D{p) by the simultaneous 
solution of eqs. and H12(l . The diffusion mean free path for other nuclei of charge Z is l{p/Z) if it is l{p) for protons. 

Let us estimate the effect of wave damping. Assuming that the cosmic-ray energy density is about 1 eV/cm'^, the 
diffusion coefficient D ~ i y. 10^* cm^/s at 1 GeV, T4 = 10 km/s, B = 5 fiG, and Cm = 1, the second term in the 
denominator in eq. (|12|l is about unity at GeV energies and falls at higher energies. We conclude that the Kraichnan- 
type cascade is significantly affected by damping on cosmic rays, and this should lead to the modification of cosmic-ray 
transport at energies less than about 10 GeV/n. 

4. SIMPLE SELFCONSISTENT MODEL 

To demonstrate the effect of wave damping, let us consider a simple case of one-dimensional diffusion with source 
distribution q = qQ{p)S{z) (corresponding to an infinitely thi n disk of cosmic-r ay sources located at the Galactic 
mid-plane z — 0) and a flat cosmic-ray halo of height H, see iJones et al.l l|2001|) . The source spectrum at i? < 40 
GV is Qo oc where 7s = 2.5 approximately. Considering energetic protons, we ignore energy losses and nuclear 

fragmentation (pioss = 0, 1/t — 0) and assume that a Galactic wind is absent {u — 0) in eq. Let us also assume 
that stochastic reacceleration does not significantly change the energies of cosmic-ray particles during the time of their 
diffusive leakage from the Galaxy and set iiT = in eq. (0. The solution of eq. Q in the Galactic disk is then 

(13, 

The escape length (the "grammage" ) , which determines the production of secondaries during the cosmic-ray leakage 
from the Galaxy, is equal to X = figvH/{2D), where is the gas surface density of the Galactic disk. 

The simultaneous solution of eqs. (|12|l and Ijl^l) allows us to find '^{p) and D{p). Introducing the function 5'o(p) = 
qo{p)H/[2DQ{p)], which is the cosmic ray spectrum for an unmodified Kraichnan spectrum, and substituting ^> from 
eq. (|13|l into eq. H12|) . one finds after two differentiations the following equation for the ratio (p ~ Dq{p) / D{p) as a 
function of a; = p^/^: 

= -M£M^(.), (14) 

with the constraint 'p{oo) = 1. Note that the same equation is valid for the ratio 5'(p)/^'o(p). 

The function '^o/x in r.h.s. of eq. (|14() can be approximated by a power law function proportional to x^^ , b — const. 
Thus for 7s = 2.5, one has & = 3 for ultrarelativistic (w = c) and b = 11/3 for nonrelativistic (w <C c ) protons 
respectively. This allows us to find the solution of eq. H14|l in an explicit form. The diffusion coefficient D{p) = 
Do{p)/if{p) is then: 

D{P) - Do{p)—jj-—jj—-^ w{p) = ,^^^^Jgpm^,{p), (15) 

r [(6 - l)/(6 - 2)J Ji/(b_2) [2w(p)J 3(6 - 2) V 

where r(2;) is the gamma function, J^{z) is the Bessel function of the first kind. With p decreasing from infinity and 
the corresponding increase of w{p) from zero, the Bessel function goes to zero at some wij)^) — (where — 1.92 at 
6 = 3, and = 1.64 at 6 = 11/3). This means that D{p) becomes infinite at p — p^ because of complete termination 
of the cascade; D{p) = oo at, p < p^,. 

The expansion of the Bessel function for small and large arguments results in the following approximations for 
eq. ((TKl) : 

^^^^ " l-{r[(6-l)/(6-2)r/rK26-3)/(6-2)]}.^(p) ^^^^ « ^''^ 

and 

D{p)^ Do{p) r'. . ^^TT at 1 < w(p) < w,. (17) 

F [(6 - l)/(6 - 2)] sm [f (36 - 4)/(6 - 2) - 2w{p)\ 

The last approximate expression for D{p) has a singularity at w{p) = = (7r/8)(36 — 4)/(6 — 2); = 1.96 at 6 = 3, 
and — 1.65 at 6 = 11/3 which is close to the values obtained from eq. H15|) . The corresponding p^ can be found 

from the equation p*^^5'o(p*) — [37r(36 — 4)/16]^5^^. 

Eqs. (|15|) - (|17|l exhibit the rapid transition from unmodified to infinite diffusion. Fig. ^ demonstrates that eq. H15|l 
derived in the considerably simplified model of cosmic-ray propagation (dash-dot) agrees well with the shape of the 
diffusion coefficient found in the numerical simulations (solid) described in the next Section. The results of this Section 
can be used for approximate estimates of cosmic ray diffusion coefficient without invoking the full-scale calculations 
based on the GALPROP code. 
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5. NUMERICAL SOLUTION WITH GALPROP CODE 

The numerical calculations have been made using the GALPROP code for cosmic-ray propagation. GALPROP is a 
flexible numerical code written in C-l — h which incorporates as much realistic astrophysical input as possible together 
with latest theoretical developments. The model is designed to perform cosmic-ray propagation calculations for nuclei 
(isotopes of H through Ni), antiprotons, electrons and positrons, and computes 7-rays and synchrotron emission in the 
same framework. The current version of the model includes substantial optimizations in comparison to the older model 
and a full 3-dimensional spatial grid for all cosmic-ray species. It explicitly solves the f ull nuclear reaction network 
on a spati ally reso lved grid. GALPROP models have been described in detail elsewhere ijStrong fc MoskalenkollTOQ^ 
iMoskalcnk o et an i2002. 2003b). For the present calculation we use three-dimensional Galactic models with cylindrical 
symmetry; (R,z,p) - spatial variables plus momentum. 

GALPROP solves the transport equation eq. Q for all cosmic-ray species using a Crank-Nicholson implicit second- 
order scheme. The spatial boundary conditions assume free particle escape. The diffusion coefficient as a function 
of momentum and the reacceleration or convection parameters are determined by boron-to-carbon (B/C) ratio data. 
The spatial diffusion coefficient is taken as I? = k(3{R/ Rq)'^ , k = const, (3 = v/c, if necessary with a break (a = oi 
below rigidity Rq, a — a2 above rigidity Rq). The injection spectrum of nucleons is assumed to be a power law in 
momentum, q{p) cx and may also have breaks. To account for stochastic reacceleration in the interstellar medium, 
the momentum-space diffusion coefficient K is related to the spati al coeHicient D via the Alfven speed Va, see eq. 0. 

The interstellar hydrogen distribution l)Moskalenko et aT]l2002D uses H i and CO surveys and information on the 
ionized component; the helium fraction of the gas is taken as 0.11 by number. The H2 and H i gas number densities in 
the Galactic plane are defined in the form of tables, which are interpolated linearly. The conversion factor is taken as 
Xco= NiijWco = 1.9 X 10^" mols. cui-'^/(K km s^^) (Strong & Mattox 1996). The extension of the gas distribution 
to an arbitrary height above the plane is made using analytical approximations. The halo size is assumed to be -ff = 4 
kpc. The distribution of cosmic-rav sources is chosen to reproduce the cosm ic-ray distribution determined by analysis 
of EGRET 7-ray data l|Strong fc Mattoxl[l996t iStrong fc MoskalenkollT99^ : 

.(A)\.,(_,,ii_Jk__j£l_), 

where qo is a normalization constant, Rq = 8.5 kpc, a = 0. 5 and = 1.0 are par ameters. We note t hat the adapted 
source distribution^ is flatter than the distribution of SNR l)Case fc Bhattac harva 199^ and pulsar s dLorimcr 2004). 

The code includes cross-section measurements and energy dependent fitting functions ({Strong fc Moskalcnko 2001). 
The nuclear reaction network is built using the Nuclear Data Sheet s. The isotopic cross section database is built 
using the extensive T16 Los Alamos compila tion of the cross sections (iMashnik et al.i ri998) and modern nuclear codes 
CEM2k and LAQGSM fMa shnik et al.ll2004|l . The most important isotopic production cross s ections (^H, '^H, ^He, Li 
Be, B, Al, CI, Sc, Ti, V, Mn) are cal culated using our fits to major production channels fe.g.. lMoskalenko et alJl2'00lL 
I2003bt IMoskalenko fc Mashnikll2003|) . Other cross sections are calculated using phcnomenological approximations by 



^ A recent analysis has shown that the apparent discrepancy between the radial gradient in the diffuse Galactic 7-ray emissivity and the 
distribution of SNR, believ ed to be the sources of cosmic rays, can be plausibly solved by adopting a conversion factor Xqo which increases 
with R (see a discussion in lStrong et al.l|2004l) . 
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Fig. 3. — B/C ratio as calculated in plain diffusion model (PD model), reacceleration model (RD model), and diffusive reaccelera- 
tion with damping model (DRD model). Lower curve - LIS, upper — modulated (<E> = 450 MV). Data below 200 MeV/nucleon; ACE 
Ipavis et al. 2000), Ulysses iDuVernois. Simpson. & Thaver 1996), Voyager {.L ukasiak. McDonald. & Webber 199^; high energy data: 
HEAO-3 tEngelmann et al.M199n.1 . for other references see [Stephens fc Streitmatten lil99!'^V 



I Webber et "al. (1990) (code WNEWTR . FOR versions of 1993 and 2003) and/or Silberberg and Tsao fSilbcrb erg et alJll99S() 
(code ^YIELDX_0 11000 .FOR version of 2000) renormalized to the data whe re it exists. For p A inelastic c ross section we 
adap ted the parametrization by Barashcnkov & PolansO (code CROSE C.lBarash enkovll993: Barashcn kov fc Polanskil 
Il994|) . The details of proton and antiproton cross sections are given in Moskalenko et al. (2002). Secondary positron 
and electron production is computed using the formalism described in Moskalenko & Strong (19.98i), which includes 
a reevaluation of the secondary tt^- and iC^-meson decay calculations. Energy losses for nucleons and leptons by 
ionization, Coulomb scattering, bremsstrahlung, inverse Compton scattering, and syn chrotron are included in the 
calculations. The heliospheric modulation is treated using the force-field approximation (|Gleeson fc AxfordllTOfil . 

The reaction network is solved starting at the heaviest nuclei (i.e., ^''Ni). The propagation equation is solved, 
computing all the resulting secondary source functions, and then proceeds to the nuclei with A — 1. The procedure 
is repeated down to A = 1. In this way all secondary, tertiary etc. reactions are automatically accounted for. To be 
completely accurate for all isotopes, e.g., for some rare cases of /3^-decay, the whole loop is repeated twice. 

For the present work, a new iterative procedure which includes a self-consistent determination of the cosmic-ray 
diffusion coefficient was developed. In the first step, cosmic-ray propagation is calculated using the undisturbed 
diffusion coefficient Do{p). In the second, we use the propagated proton spectrum at every spatial grid point to re- 
calculate the modified diffusion coefficient according to eq. p2(l . In the third step, cosmic-ray propagation is calculated 
using the new diffusion coefficient. Steps 2 and 3 are repeated until convergence is obtained. The convergence is fast 
and the procedure requires only few iterations. In this way, using the secondary to primary ratio B/C, one can derive 
the damping constant g = 0.085 which is considered as an adjustable parameter. 

The new self-consistent model of diffusive reacceleration with damping (DRD model) is compared with two reference 
models which were earlier found to give good fits to cosmic ray data: the plain diffusion model (PD model) with an 
ad-hoc break in the diff usion coefficient and the model with distributed reacceleration (DR model) and power-law 
diffusion with no breaks ijMoskalenko et aDl2002D . The parameters of these models found from the fit to the observed 
spectra of primary and secondary nuclei are summarized in Tabled Note that the parameters of the DRD model given 
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TABLE 1 
Propagation parameter sets. 





Injection index'' 


Break 


Diffusion coefficient @ 3 GV 


Alfven speed. 




Model 


Nucleons, 7^ 


Electrons, 7^ 


rigidity, GV 


K, cm^ s ^ 


Index, a 


Va, km s-i 


galdef-file 


Plain Diffusion (PD) 


2.30/2.15 


2.40 


40 


2.2 X 10^8 


0.0/0.60*= 




44.999726 


Diffusive 

Reacceleration (DR) 


1.80/2.40 


1.60/2.50 


4 


5.2 X 1028 


0.34 


36 


44.599278 


Diffusive 

Reacceleration with 
Damping (DRD) 


2.40/2.24 


2.70 


40 


2.9 X 10^8 


0.50 


22 


44.999714kr 



Note. — Adopted halo size H — 4 kpe. 
^Index below/above the break rigidity. 
''Index below/above Rq = 3 GV; D = /3"^K(i?,/flo)°- 




10"'' 10° 10^ 10^ 10^ 
Kinetic energy, GeV/nucleon 

Fig. 4. — Spectrum of carbon calculated in plain diffusion model (PD model), reacceleration model (RD model), and diffusive reacceleration 
with damping model (D RD model). Upp er curves — LIS, lo wer curves modulated using fo rce field approximation ($ = 55 MV). Data: 
ACE IDavis et ani2000l IZDOlD . HEAO-3 lEngelmann et al.lll990D . for other references see IStephens fc Streitmattej 119981) (symbols are 
changed) . 

in Tableware the "seed" values used in the iteration procedure. The diffusion coefficient after iteration calculations 
is shown in Fig. 

The diffusion coefficients in all three models are presented in Fig. ^ The DR model diffusion coefficient has a weak 
energy dependence with a single index 1/3 as dictated by the assumed Kolmogorov spectrum of interstellar turbulence. 
The PD model requires a break in the diffusion coefficient and an additional factor of to be able to match the B/C 
ratio at low energies. The DRD model diffusion coefficient lies between these at high energies and has a sharp increase 
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Fig. 5. — Proton spectra as calculated in plain diffusion model (PD model), reacceleration model (RD model), and diffusive reacceleration 
with damping model (DRD model). Upper curve - LIS, lower - modulated to 550 MV. Thin dotted line shows the LIS spectrum best 
fitt ed to the data abov e 20 GeV (Moskalenko et al. 2002). Data: AMS (Alcaraz et al. 2000b), BESS 98 iSanuki et al. 2000), CAPRICE 
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at 1 GV. Of course, the actual mean free path length can not be infinite. We put an upper limit of ^15 pc based on the 
estimate of streaming instability effect that arises at low rigidities and leads to the generation of additional turbulence 
which limits cosmic ray escape from the Galaxy, see discussion below in Section 6. The selfconsistent spectrum of 
waves calculated in the DRD model is shown in Fig.|21 The spectrum bends downward at wavenumbers fc > 6 x 10^^^ 
cm~^ due to the wave damping on cosmic rays. 

The B/C ratio and Carbon spectra (primary) after modulation look almost identical in all three models (Figs.|3121). 
The modulated proton spectrum (Fig. matches the data in all three cases after some tuning, while the interstellar 
spectra are different. The most dramatic difference is the reduction of proton flux in the DRD model at low energies, 
consistent with the steep increase of the diffusion coefficient. The antiproton spectrum (Fig. |HJl appears to be very 
sensitive to the model assumptions and might help to discriminate between the models. The DR model produces 
too few antiprotons, a well known effect (Moskalenko ct al. 2002). The PD and DRD models are both consistent 
with antiproton measurements in the heliosphere, but predict very different spectra in the interstellar medium; the 
interstellar antiproton flux at low energies (<600 MeV) in DRD model is an order of magnitude lower than in two other 
models. In both reacceleration models (DR and DRD), the majority of low-energy antiprotons come from inelastic 
scattering (so-called "tertiary" antiprotons). 

Figs. [3 and |S1 show secondary positrons and primary plus secondary electrons as calculated in all three models. The 
spectra are similar in the PD and DR models, while DRD spectra exhibit lower intensities at low energies. This may 
be an observable effect since the models predict different synchrotron emission spectra (electrons). 

6. DISCUSSION 

Damping on cosmic rays may terminate the slow Kraichnan-type cascade in the interstellar medium at fc lO^^'^ 
cm^^. Our estimates were made for the level of MHD turbulence which produces the empirical value of cosmic-ray 
diffusion coefficient. This finding suggests a possible explanation for the peaks in secondary/primary nuclei ratios 
at about 1 GeV/n observed in cosmic rays: the amplitude of short waves is small because of damping and thus the 
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low energy particles rapidly exit the Galaxy without producing many secondaries. There is no other obvious reasons 
for a sharp cut off in the wave spectrum. If the concept of MHD turbulence by iGoldreich fc Sridhad l|1995D works 
for interstellar turbulence, the MHD waves we are dealing with in this context are the fast magnetosonic waves. The 
Alfven waves propagate predominantly perpendicular to the magnetic field and because of this they do not significantly 
scatter cosmic rays. It also explains why radio scintillation observations show no sign of the termination of el ectron 
density fluctuations at wave numbers between 10^^^ to 10~^ cm~^. According to lLithwick fc GoldreichI 1)200 Ij) these 
fluctuations are produced by the slow magnetosonic waves with k± ^ /cj| which are almost not damped on cosmic 
rays. An alternative explanation is that the wave damping on cosmic rays and the radio scintillations mainly occur 
in separate regions of the interstellar medium (see below in the discussion on the "sandwich" model of cosmic ray 
propagation in the Galaxy) . Some minor contribution to the observed scintillations is possible from fast magnetosonic 
waves interacting with energetic particles, and in this respect it is of interest that the observations may need an 
enhancement in the power on large "refra ctive" scales 10^^ — lO^'' cm relative to the power on small "diffractive" scales 
10^ — 10^" cm (jLambert fc Rickettll2000l) . This may indicate the cutoff of the spectrum of fast magnetosonic waves 
due to cosmic ray action. 

While the mere fact of a wave spectrum steepening under the action of damping on cosmic rays can be described 
by a simple eq. © with some characteristic time for nonlinear wave interactions Tni[W{k), k], the exact form of the 
function W{k) at large k where the damping is significant depends critically on the form of the equation for waves. It 
involves in particular to the vanishing of W{k) at some fc* (~10~^^ cm~^) and the corresponding singularity of D(j)) 
at some p» (~1 GV) found in our calculations. Less significant in this sense is our approximation of the resonant wave 
number in eq. (jSJ, which does not include the particle pitch angle in an explicit form. (If it were included, the term 
with g in eq. [5| would have a third integration over the pitch angle.) We note however that this effect was included 
in the derivation of eq. ^Hl for the attenuation rate. 

In the context of the approximations adopted in the present work, the problem of cosmic-ray transport at p < 
arises. The free streaming of cosmic rays from the Galaxy leads to an instab ility and to the growth of waves which 
scatter particles and thus slow down the streaming, see e.g. iBerezinskii et al.] ()199niD . We shall consider the processes 
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Fig. 7. — Positron flux as calculated in plain diffusion model (PD model), reaccelerati on model (RD model) , and diffusive reacceleration 
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l2000h . HEAT 94-95 IDuVernois et aljl200ll) . MASS 91 IGrimani et'ani200l) . 



at low energies in a separate work. Here we note only that the streaming instability develops above the Galactic disk 
at -D(l GV) > Vs.il ~ 3 X 10^^ cm^/s (if the wave damping is absent) and leads to the difFusive-convective transport of 
cosmic rays. The given estimate of the diffusion coefficient at 1 GV follows from the condition of cosmic ray streaming 
instability Ucr > Va where Ucr is the bulk velocity of cosmic ray gas. The bulk velocity is Ucr ~ ScrC where the cosmic 
ray anisotropy perpendicular to galactic disk is Scr ~ D/cH. The Alfven velocity in galactic halo is about 2 x 10^ 
cm/s. It is important that the magnetic rigidity 1 GV corresponds to a kinetic energy 0.43 GeV for protons, and 0.13 
GeV/n for nuclei with charge to mass ratio Z/A = 1/2. The Galactic spectrum of cosmic-ray protons and nuclei at 
such low energies can not be derived from direct observations at the Earth because of strong modulation in the solar 
wind. 

Dissipation other than on cosmic rays has been neglected in the present work, though it may completely destroy the 
MHD casca de or con siderably change its angul ar dis tribution and thus affect D(p) in a large part of the interstellar 
medium, see iMcIvorl fe77), Cesarskv (1980), lYan fc LazarianI l)2004l) . The region of the cosmic ray halo is the most 
"safe" in this sense l|Yan fc Lazariami2004) . In particular, dissipation on ion- neutral collisions may destroy MHD 
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turbulence in the Galactic disk but not in the halo where neutrals are absent. We then come to the "sandwich" model 
for cosmic-ray propagation, with different diffusion coefficients in the Galactic disk, Dg, and in the halo, Dh (it is 
assumed here that some scattering is present in the Galactic disk and the diffusion approximation works). In this 
case, the mean matter thickness traversed by cosmic rays does not depend on diffusion inside the disk even at zero 

gas density in the halo: X = ^ (^2S~ ~ tJ-g^H / (2D\y) ijGinzburg fc Ptuskinlll976fl . The energy dependence 

of secondary to primary ratios in cosmic rays is determined by the diffusion coefScient in the cosmic ray halo where 
the model developed in the present paper is applied. 

The estimate based on the empirical value of the diffusion coefficient for GeV particles (see Section 1) gives the level 
of turbulence at the principal scale SB^^^/B^ ~ 0.03 for a Kraichnan-type spectrum W{k) oc k^^^^, and SBf^^/B^ ~ 1 
for a Kolmogorov-type spectrum W{k) cc k'^^^, if ki^ = 10~^^ cm~^. At the same time, the data on Faraday rotation 
angles favor the Kolmogorov spectrum with 6B^^^/B^ ~ 1 and = 10~^^ cm~^. The cascades of Alfven waves (with 
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the scaling fc~^/^) and the fast magnetosonic waves (fc^^^^) are independent in the lGoldreich fc Sridhad l)1995() model 
of MHD turbulence, and the amplitude of Alfven wave cascade may dominate at the principle scale. Also, in the 
"sandwich" model described in a previous paragraph, the turbulence, which determines the confinement of cosmic rays 
in the Galaxy, is distributed in the halo of size H ^ 4 kpc whereas the observations of interstellar turbulence refer to 
the Galactic disk and the adjacent region, where it can be much stronger. 

7. CONCLUSIONS 

On the whole, the empirical diffusion model for cosmic rays with energies from 10^ to 10^^ eV implies the presence of 
random magnetic field with an extended power law spectrum of fluctuations W{k) cc A:"^"*"", a < 0.5 at wave numbers 
from 3 X 10~^^ to 10"^" cm~^. The existence of such a turbulence spectrum in the interstellar medium seems confirmed 
by various astronomical observations. It should be emphasized that this does not prove the existence of a spectral 
cascade or spectral transfer throughout this enormous wavenumber range, however tempting that conclusion may be. 
Among other puzzles, an oddity is the absence of a spectral feature on spatial scales where ion-neutral coUisional 
processes should be most pronounced. The two special cases of turbulence spectrum with a = 1/3 and a — 1/2 which 
correspond to the Kolmogorov and the Kraichnan spectra respectively are used in popular versions of the diffusion 
model as described in Section 5. In both cases, as the present work has demonstrated, one can get a satisfactory fit 
to the data on energy spectra of secondary and primary nuclei if account is taken for wave damping on cosmic rays 
for a slow Kraichnan cascade. The fit can be obtained either in the DR model with reacceleration on the Kolmogorov 
spectrum with no significant effect of cosmic ray damping (because the Kolmogorov cascade is fast) or in the DRD 
model with relatively weak reacceleration on the Kraichnan spectrum which is significantly modified by cosmic ray 
damping. Some problems still remain to be solved. The main difficulties with the Kolmogorov spectrum are, firstly, 
the contradiction with the leading theory of MHD turbulence where this spectrum is associated with perpendicular 
propagating Alfven waves which almost do not scatter cosmic-ray particles; and secondly, the low flux of antiprotons 
characteristic of models with relatively strong reacceleration. A major problem of concern for diffusion on a Kraichnan 
spectrum is the relatively strong dependence of diffusion on energy which leads to an unacceptably large anisotropy 
of cosmic rays especially above lO^"* eV, seeiJ.oncs ct al. (2001), Ptuskin et al. (_2QQ.3ii)- 

Let us emphasize again that the models of cosmic ray propagation discussed in the present paper assume that 
the MHD turbulence re quired for cosmic-ra y scattering is pro duced by some external sources. An alternative is the 
Galactic wind model bv iZirakashvili et al.l (^96) and Ftuskir Tet al.l (^97). In this model, the cosmic-ray pressure 
drives a wind with a frozen-in regular magnetic held which is shaped into huge spirals (the radius of the Galactic wind 
cavity is about 300 Kpc). The MHD turbulence in the wind is created by the streaming instability of cosmic rays 
moving predominantly along the regular magnetic held lines outward from the Galactic disk. The level of turbulence 
is regulated by the nonlinear Landau damping on thermal ions. The model explains well the cosmic-ray data up 
to ultra-high energies ~10^^ eV with the exception of the observed low anisotropy (about the same difficulty with 
anisotropy as occurs in the diffusion model with a given Kraichnan spectrum). As regards interpretation of cosmic-ray 
observations, a preference cannot yet be given to any of the models of cosmic-ray transport in the Galaxy discussed 
above. 
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